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(N 

t-H We present 2+1 flavor Lattice QCD calculations of the nucleon scalar and tensor charges. Using 

the BMW clover-improved Wilson action with pion masses between 149 and 356 MeV and three 
£Nj source-sink separations between 0.9 and 1.4 fm, we achieve good control over excited-state contam- 

ination and extrapolation to the physical pion mass. As a consistency check, we also present results 
O from calculations using unitary domain wall fermions with pion masses between 297 and 403 MeV, 

and using domain wall valence quarks and staggered sea quarks with pion masses between 293 and 

597 MeV. 

ON 

t— I PACS numbers: 12.38.Gc, 12.60.-i, 13.30.Ce, 14.20.Dh 

I. INTRODUCTION 

p , 

A quantitative understanding of the quark structure of hadrons is essential to relate experimentally observed 
hadronic properties to fundamental processes occurring at the quark level. Thus, Lattice QCD can play a crucial 
role in determining parameters of the Standard Model from experiment and in predicting the experimental effects of 
interactions beyond the Standard Model. 

One example is the search for the effects of new scalar and tensor couplings beyond the familiar weak interactions 
of the Standard Model in the decay of ultra-cold neutrons as discussed in detail in pp. These couplings would 
contribute through the matrix elements (p\ud\n) and {p\ua ^ v d\n) . To leading order in the recoil approximation, these 
matrix elements are proportional to the isovector scalar and tensor charges gs and gx of the nucleon. Hence, Lattice 
QCD predictions of these charges are crucial for determining the sensitivity of the decay process to the non-standard 
\£) couplings. 

Another example concerns the search for an electric dipole moment of the neutron, which constrains the possible 
values of the CP-violating strong interaction 9 angle. The constant of proportionality between the electric dipole 
moment and 9 can be calculated directly using a background electric field [5J[3] or from the CP-odd form factor F3, 
J> cf . [3J [5] . However, it is also instructive to relate the neutron electric dipole moment near the chiral limit at small 9 
to the CP- violating pion-nucleon coupling constant g-xNN, cf. [BJ. This assumes that |-7riV) intermediate states provide 
the dominant contribution to the dipole moment, which is the case as m n — > 0. In turn, g-^NN is proportional to 
the isovector scalar charge gs of the nucleon [BJ, which thus again appears as an important parameter influencing 
the search for new physics. Given that the physical limit occurs at finite m n , where there will be corrections to the 
relation between the electric dipole moment of the neutron and g^NN, a comparison of this approximation to direct 
evaluations can also provide insight into the importance of these corrections. Another important effect of g^NN is the 
fact that it induces CP violation in nuclear forces resulting in strongly enhanced nuclear electric dipole moments in 
certain cases [7j. 

These two examples directly motivate the present work in which we calculate the nucleon's isovector charges gs 
and gx- Since the computationally formidable disconnected contributions from couplings to the sea quarks cancel in 
these isovector charges, we only calculate the connected contributions. Extensive Lattice QCD data generated with 
three different lattice actions at pion masses ranging between 149 MeV and 597 MeV are reported, together with their 
extrapolations to the physical pion mass. 

In addition to the primary motivation described above, our results also contribute to a more detailed understanding 
of the nucleon sigma term cr/v = ^(m u + md)(N\uu + dd\N) , which determines the part of the nucleon mass generated 
by the light quark degrees of freedom via the spontaneous breaking of chiral symmetry. Whereas we do not have 
the resources to evaluate the disconnected contributions to <7n at present, the connected contributions to (N\uu\N) 
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Label 


a [fm] 


L x x L t 


m n [MeV] m 


7T L x a 


# confs # 


meas 










Wilson-clover 






Wl 


0.09 


32 3 


x 64 


317(2) 


4.6 


103 


824 


W2 


0.116 


48 3 


x 48 


149(1) 


4.2 


646 


7752 


W3 


0.116 


32 3 


x 48 


202(1) 


3.8 


457 


5484 


\H T A 

W4 


0.116 


32 3 


x 96 


253(1) 


4.8 


202 


2424 


W5 


0.116 


32 3 


x 48 


254(1) 


4.8 


420 


5040 


W6 


0.116 


24 3 


x 48 


254(1) 


3.6 


418 


10032 


W7 


0.116 


24 3 


x 48 


303(2) 


4.3 


128 


768 


W8 


0.116 


24 3 


x 48 


356(2) 


5.0 


127 


762 










Domain wall 








Dl 


0.084 


32 3 


x 64 


297(5) 


4.0 


615 


4920 


D2 


0.084 


32 3 


x 64 


355(6) 


4.8 


882 


7056 


D3 


0.084 


32 3 


x 64 


403(7) 


5.5 


527 


4216 


D4 


0.114 


24 3 


x 64 


329(5) 


4.6 


399 


3192 










Mixed action 








A T 1 
Ml 


n 1 1 a 
U.124 


20 3 


x 64 


293(6) 


3.7 


a fin 


3680 


M2 


0.124 


28 3 


x 64 


356(7) 


6.3 


272 


2176 


M3 


0.124 


20 3 


x 64 


356(7) 


4.5 


628 


5024 


M4 


0.124 


20 3 


x 64 


495(10) 


6.2 


483 


3864 


M5 


0.124 


20 3 


x 64 


597(12) 


7.5 


562 


4496 



Table I. Lattice ensembles used for scalar and tensor charge calculations. 



and (N\dd\N) can be given separately and are expected to furnish most of the strength of the nucleon sigma term; 
cf. [5HTU] for recent lattice studies focusing on ajy. These quantities are also relevant to dark matter searches based 
on Higgs-mediated couplings of baryonic to dark matter [TTrfTB"] . For example, the neutralino-nucleon scalar cross- 
section considered in [TJJ can be cast as a linear combination of terms proportional to (N\uu\N) and (N\dd\N). Thus, 
again, the scalar charge gs of the nucleon (including both isoscalar and isovector components) is useful in quantifying 
searches for these dark matter candidates. 



II. LATTICE METHODOLOGY 



The full set of lattice ensembles on which we have calculated our observables is listed in Tab. |U We make use of 
three different lattice actions. 

Our main results are calculated on ensembles at light pion masses, with tree-level clover-improved Wilson fermions 
coupled to double HEX-smeared gauge fields, as used by the BMW collaboration [T3]. Operators are renormalized 
nonperturbatively using the Rome-Southampton method. For the scalar charge, we use renormalization factors cal- 
culated by the BMW collaboration [T3] , and for the tensor charge we performed our own calculation, see Appendix [A| 
These factors are given in Tab. [TTJ 

a [fm] Z s Z T 



0.09 1.107(16)(22) 1.011(5) 
0.116 1.115(17)(30) 0.9624(62) 



Table II. MS renormalization factors at /x = 2 GeV for the tensor and scalar charge on the Wilson-clover ensembles. 

For comparison and for a consistency check, we also present results on ensembles with two different lattice actions 
used for earlier nucleon structure calculations. First, we use unitary domain wall quarks on ensembles generated 
by the RBC and UKQCD collaborations [IB] US] ■ Details of our analysis methods are in Ref. [T7] . Operators are 
renormalized nonperturbatively using the Rome-Southampton method. Zg 16J and, on the coarse ensemble, Zt [TBI 

were calculated by the RBC collaboration. On the fine ensembles, we found Z^ IS< " 2 GeV ' = 0.8168(9), see Appendix A 
Second, our mixed-action scheme [191 120] uses domain wall valence quarks on gauge configurations with Asqtac 
staggered sea quarks generated by the MILC collaboration [21]. Renormalization factors are calculated in the same 
way |22j as those used in previous calculations of nucleon generalized form factors. One-loop perturbation theory is 
used to calculate the renormalization factors for all quark bilinear operators evaluated at the one-loop coupling constant 
and /i 2 = 1/a 2 . Because HYP smearing suppresses loop integrals, the ratio •Zo,pert/-Z4,pert) where O denotes a general 
bilinear and A is the axial current, is within a few percent of unity indicating perturbative corrections are already 
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small at the one-loop level for this ratio. The multiplicative wave function renormalization in the fifth dimension 
appearing in all renormalization factors is included non-perturbatively by using the five-dimensional conserved axial 
current for domain wall fermions to calculate ^A.nonpcrt and calculating Zq = (i>e>,pert/i>A,pert) ■ Zyi,nonpert- The 
relevant matrix elements are [55] £J>i,pert = 0.964, Zs, pe rt — 0.971, and import = 0.987, and the result is evolved to 
[i = 2 GeV. 

For the scalar charge, we can alternatively make use of the renormalization group invariant combination (m s — 
m u d)gs- Dividing by the physical m s — m u d, as measured in 2+1 flavor Lattice QCD calculations in the MS scheme 
at 2 GeV [THl I2"3"H2"T] , and multiplying by a ratio of differences of pseudoscalar meson masses to cancel the leading 
dependence on quark masses, we get 



, „ , /^barc ™ bare \ ~bare rn^ rrp 1 

MS(2 GoV) _ [m s - m ud )g s m if,phys "Vphys m 
9S ~ phys phys m 2 _ ™2 ' ^ ' 

m% - m„/ m K 



where the physical rriK and in the isospin limit are from Ref. [23] and their values on our mixed-action ensembles 
were computed in an earlier work [28]. On the mixed-action ensembles, this approach yields values of gs that differ 
from the perturbatively renormalized values by between 0.4% and 1.3%. Assuming the perturbative renormalization 
of the tensor charge has similar errors, we use our perturbative approach for both gs and gx and conservatively 
estimate a systematic error due to renormalization of 2%. 

We compute nucleon forward matrix elements using the usual ratio-plateau method jTTl QjJ] . Beginning with two- 
point and three-point functions, 

C 2pt (t,P) = (N(p = P,t)N(x = 0,0)), (2) 
C 3 ° pt (r, T; P) = (N(p = P, T)(D{p = 0, t)N(x = 0, 0)) , (3) 

where N is our lattice nucleon interpolating operator, we compute their ratio, 

RO{ T p )= (4) 
C 2pt (T,P) 

At sufficiently large r and T ~ r, contributions from excited states are negligible, and the ratio gives us the matrix 
element (N(P)\0\N(P)) . In practice, for a fixed source-sink separation T, we take the average over the central two 
or three points of the plateau as the matrix element. The matrix elements give us the isovector scalar and tensor 
charges: 

(N(P)\uu - dd\N(P)) = g s u{P)u(P) (5) 
(N(P)\ua^u - da^dlNiP)) = g T u(P)a fiU u{P). (6) 

We take the weighted average of results measured with P = and P = ^(-1,0,0). The source-sink separations T 
that we use are listed in Tab. |III| On the Wilson action ensembles, we perform measurements at three source-sink 
separations in order to better identify systematic errors from excited-state contamination. However, our main results 
are presented using the middle separation, with T m 1.16 fm. 

Action a [fm] T/a 

Mixed 0124^ 9 

Domain wall 0.084 12 
Domain wall 0.114 9 
Wilson-clover 0.09 10, 13, 16 
Wilson-clover 0.116 8, 10, 12 

Table III. Source-sink separations T. 
Chiral perturbation theory results for gA and gx are summarized in Appendix [B] 



III. RESULTS 



The dependence on source-sink separation for the Wilson-clover data is shown in Fig. [T] for the scalar charge and 
in Fig. [2] for the tensor charge. Neither observable shows a strong statistically significant dependence on source-sink 
separation. 
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Figure 1. Scalar charge measured at three source-sink separations on the Wilson-clover ensembles, as enumerated in Tab.|I] 
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Figure 2. Tensor charge measured at three source-sink separations on the Wilson-clover ensembles, as enumerated in Tab. [T] 



The main results for the scalar charge are presented in Fig. [3] There is broad agreement between the different 
actions, although they are not entirely consistent. In particular, the two mixed-action ensembles at = 356 MeV 
have points that lie significantly below the point for the domain wall ensemble at = 355 MeV. We estimate the 
value of the scalar charge at the physical pion mass, 135 MeV, in three ways: 

1. Taking the value from our ensemble with the lowest pion mass, 149 MeV. Including the error in Zg, this is 
55 = 1-01(27). 
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Figure 3. gs versus raj. Error bars are purely statistical and do not include errors in renormalization factors, which are 
correlated across ensembles with the same action. Two chiral fits are shown: a three-parameter fit to the coarse Wilson-clover 
ensembles (orange) and a four-parameter fit to all shown ensembles (magenta). 



2. Performing a 3-parameter chiral fit to the set of seven a = 0.116 Wilson-clover ensembles. This is a good fit 
with x 2 /dof = 3.99/4 and is shown in Fig. [3} Extrapolating to the physical pion mass yields gs = 1.08(28). 

3. Performing a 4-parameter chiral fit to the full set of ensembles. The tension at higher pion masses between 
the domain wall and the mixed action ensembles is reflected in % 2 /dof = 30.23/13. Extrapolation yields 
.95 = 1-08(23). 

Our tensor charge results are shown in Fig. [4] Note that the vertical scale is expanded and the error bars are much 
smaller. Again there is general agreement between different lattice actions. The figure shows some inconsistency 
between the domain wall and mixed-action ensembles at larger pion masses. However, this is close to the 2% correlated 
error that we assign to the mixed-action results to account for perturbative renormalization. As with the scalar charge, 
we estimate the physical value of gx in three ways: 

1. Using the value from the m T = 149 MeV ensemble yields gx = 1.049(23). 

2. Performing a 2-parameter chiral fit to the coarse Wilson-clover ensembles. As with the scalar charge, the fit 
works well with % 2 /dof = 4.49/5. However, the fit chooses /1 2 = 0, which removes the chiral log. At the physical 
pion mass, this gives gx — 1.038(11). 

3. Performing a 3-parameter fit to the full set of ensembles. Because of the correlated error for the mixed-action 
ensembles, this fit also works well, giving % 2 /dof = 19.69/14. Extrapolating to the physical pion mass yields 
g T = 1.037(20). 



IV. SUMMARY AND CONCLUSIONS 



Motivated by the facts that contributions of non-standard scalar and tensor couplings to neutron decay are pro- 
portional to isovector gs and gx, and that g^NN, which produces parity- violating components of nucleon-nucleon 
interactions, is proportional to isovector gs, we have undertaken a comprehensive set of calculations of gs and gx 
using Lattice QCD. 

Our calculation includes a number of significant advances. We have utilized three sets of ensembles that cover 
different pion mass regions for which we use the same calculational methodology. The primary data for this analysis 
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Figure 4. gr versus m^. Error bars are purely statistical and do not include errors in renormalization factors, which are 
correlated across ensembles with the same action. Two chiral fits are shown: a two-parameter fit to the coarse Wilson-clover 
ensembles (orange) and a three-parameter fit to all shown ensembles (magenta). 



were obtained using the smeared Wilson-clover action developed by the BMW collaboration and include pion masses 
in the vicinity of 150, 200, 250, 300, and 350 MeV. The 149 MeV pion mass plays a crucial role in our analysis. It is 
so close to the physical pion mass that it virtually eliminates uncertainty from extrapolation to the physical mass and 
is well below the lowest mass used in other calculations. Furthermore, given the known strong source-sink separation 
dependence in some observables such as the momentum fraction, our careful study of 3 separations in Figs. [T] and 
Fig. [2]from 0.93 fm to 1.39 fm is important in ruling out contributions of excited states near the physical pion mass. 

To control systematic effects from finite lattice volume, we have satisfied the standard criterion m^L x a > 4 for 
nearly all of our ensembles, including the crucial m T = 149 MeV ensemble where m v L x a — 4.2, cf. Table|l] In addition, 
we carried out an explicit test for finite volume effects using the two ensembles W5 and W6 at — 254 MeV, which 
are identical except for their lattice volume. Within the statistical accuracy of our calculation, no significant deviations 
between the results obtained in these two ensembles were found. 

Comparison of values of g$ and gr for the ensembles Wl with lattice spacing a = 0.09 fm and W7 with a = 0.116 fm 
at similar values of m w and L x a, cf. Tabic |TJ indicates the absence of a-dependence within our statistical errors. 

Supplementing the results obtained at the pion mass — 149 MeV, we performed 3 and 2-parameter chiral fits to 
gs and gx for the low mass Wilson ensembles, where we expect chiral perturbation theory to be valid. These fits serve 
as a check on the consistency of the 149 MeV data point and have the potential to reduce the statistical error. The 
central values of the 149 MeV data points and fits were essentially identical on the scale of the errors, and it turned 
out that the error in the scalar charge was essentially unchanged while the error on the tensor charge was cut in half. 
We regard the chiral fits to the Wilson-clover data as our definitive results for the central values and statistical error. 
We also estimate the systematic error from excited-state contamination as the average of two absolute differences; 
those between the central value from the middle source-sink separation on the 149 MeV ensemble and the central 
values from the other two source-sink separations on that ensemble, i.e., \{\g§T Wa — 9x =Sa \ + \dx =10a ~ 9x =12a \)- 
This yields g s = 1.08 ± 0.28 (stat) ± 0.16 (syst) and g T = 1.038 ± 0.011 (stat) ± 0.012 (syst). We emphasize that our 
results do not hinge decisively on chiral perturbation theory. The chiral fits to the data at higher pion masses rather 
serve to buttress the calculations at m n = 149 MeV, which by themselves represent determinations of gs and gx with 
essentially no residual uncertainty due to chiral extrapolation. 

The range of pion masses for the Wilson action overlaps the range 297 to 403 MeV for the domain wall action and 
the range 293 to 597 MeV for the mixed action, so simultaneous analysis provides a valuable consistency check of 
the actions and normalization. In the region of overlap between the Wilson-clover and domain wall calculations, the 
results agreed well, confirming the underlying consistency of these two calculations including normalization. There 
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is a certain tension between the mixed action and domain wall results at higher pion masses, the origin of which we 
have not been able to pinpoint; an underestimate of the renormalization uncertainties in the mixed action case is one 
possible source. When the error in the normalization is taken into account, the 4 and 3-parameter chiral fits to g$ 
and <?t for all three sets of ensembles are consistent with the chiral fits to the low mass Wilson ensembles, as reflected 
in the close agreement of the error envelopes at the physical pion mass. 

Concerning other systematic uncertainties, we have yet to perform extensive low mass Wilson calculations at smaller 
lattice spacings to extrapolate to the continuum limit and at larger volumes to extrapolate to the infinite volume limit. 
At the present level of statistics, our comparison of two different lattice volumes at the pion mass m n = 254 MeV and 
two different lattice spacings at ~ 300 MeV yielded no evidence of significant volume or lattice spacing effects. 
However, in view of past experience with the well-studied case of the nucleon axial charge gA, future more statistically 
accurate calculations may yet necessitate more detailed investigations of a range of volumes and lattice spacings. 
Comparable calculations of gA by several groups including our own lie below the experimental value by more than 
the statistical errors, and it is likewise possible that the aforementioned or other systematic errors in gs and <?t are 
still appreciable. 
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Appendix A: Non-perturbative renormalization of the tensor charge operator 
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Figure 5. Extrapolations of the scale-independent (SI) ratio (Zf. 1 //?^*) for coarse clover-improved Wilson fermion lattice 
ensembles (W3, W5) (left) and fine domain wall fermion lattices (Dl, D2, D3) (right). 



In this section, we briefly describe the method and present the results of computing the nonperturbative renormal- 



MS(2 GeV) 



ization factors Z, 
domain wall fermion ensembles. 



for the tensor charge gr, Z T 



MS(2 GeV) lat 
i/T 



^MS(2 GeV)^ clover-improved Wilson and 



Following Ref. 



we compute amputated lattice Green functions LT^ of tensor charge operator and plane wave 
= /i 2 in the Landau gauge. We then extract the lattice renormalization factor defined 
where n tree is the corresponding tree-level amputated Green function and Z q is the quark 
field renormalization. In order to match lattice-regulated and perturbatively renormalized operators, we compute 
scale-independent (SI) renormalization factors Zf. 1 



quarks q(x) 
as {Z/Z q ) ut Ii 



- e lpx , v = 

flat jjtrec 



vSI 

Zjrp 



lat / 



Z RI'(2 GeV) v , ( ]( 
yRI'(p) 



m q — >-0, (J.— >0 



(Al) 



where the perturbative operator renormalization is calculated to order 0(otg) |31j and Z^ is the renormalization of 
the quark charge operator on the lattice, 



(Z v /Z q ) ut ■ (N(P)\v?u - Jd\N(P)) Ut = 1 . 



(A2) 



We extrapolate linearly first in m q — ¥ 0, and then in /i — !• in the range /i 2 « 6 ... 15 GeV , see Fig. [5] Finally, we 
obtain the renormalization factors for converting lattice results into the MS(2 GeV) scheme: 



,MS(2 GeV) _ / Z T 



z 



lat 



Z MS{2 GeV) ^ pcrt 
r 7 RI'{2 GeV) 



(A3) 



Appendix B: Chiral Perturbation Theory 



For chiral extrapolation of the scalar charge, we make use of the Feynman-Hellmann theorem to obtain 



9s = (p\uu - dd\p) = 



d d 
dm u duid 



(Bl) 



Strong isospin splitting of nucleon masses has been computed in SU (2) heavy baryon chiral perturbation theory, with 
explicit A degrees of freedom, to next-to-next-to-leading-order, i.e. to 0(m q ) [32]. This leads to the expression 



9 s (ml) = 2a M 



2(4^) 2 



ml log ( 2a M (6<?i + 1) + J(m n ,A, fJ,)8g 2 AN (a M 



8(g 2 A a M + 9 2 AN (a M + §7m)) - K + K ) 



(B2) 



j7mJ 



where 



7/i I 

J(m, A, /i) = (m 2 - 2A 2 ) log + 2AV A 2 - m 2 log 



A- x/A 2 " 



A + s/E^ 



2A 2 log 



4A 2 
7^" 



(B3) 



Here F^, and A are the pion decay constant, the nucleon axial charge, and the A(1232)-nucleon mass splitting, 
respectively, in the SU(2) chiral limit. The convention used here is that the physical value of F„ is approximately 
92 MeV. We fix these three parameters to values used for chiral extrapolation in a previous nucleon structure calcu- 
lation UHL 



86.2 MeV, 



9A = 1-2, 



0.2711 GeV, 



(B4) 



leaving 3 independent parameters to which we fit our 55 (m 2 ) data. If we fix only A, then the resulting fit has 4 
independent parameters. 

Chiral extrapolation formulas for the tensor charge are given in [33] . Full formulas including A loops are given via 
a number of integrals, and then it is shown that these are approximated well by an expression that includes only the 
leading non-analytic term. Also including terms that connect to the heavy quark limit: 



griml) = 6a ( 1 + <5c LNA m 2 log 



5b- 



(B5) 
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Scalar 


Tensor 


Ensemble 


m w [MeV] 


u — d 


u + d 


u — d 


u + d 


Wl 


317(2) 


1.017(248) 


5.618(388) 


1.047(38) 


0.700(46) 


W2 


149(1) 


1.006(271) 


10.127(1.501) 


1.049(22) 


0.652(33) 


W3 


202(1) 


1.009(235) 


6.632(620) 


1.052(30) 


0.622(46) 


W4 


253(1) 


1.335(262) 


7.808(906) 


1.014(20) 


0.598(55) 


W5 


254(1) 


0.901(110) 


6.027(251) 


1.041(16) 


0.676(22) 


W6 


254(1) 


0.757(195) 


5.241(935) 


1.068(28) 


0.668(33) 


W7 


303(2) 


0.622(260) 


4.923(395) 


0.995(40) 


0.656(63) 


W8 


356(2) 


0.817(168) 


5.128(268) 


1.051(35) 


0.691(46) 


Dl 


297(5) 


0.964(85) 


8.933(195) 


1.050(12) 


0.779(19) 


D2 


355(6) 


0.854(43) 


8.114(102) 


1.080(8) 


0.807(11) 


D3 


403(7) 


0.941(35) 


7.903(81) 


1.087(8) 


0.802(11) 


D4 


329(5) 


0.839(65) 


8.193(153) 


1.087(12) 


0.827(18) 


Ml 


293(6) 


1.012(208) 


6.624(278) 


1.069(27) 


0.543(34) 


M2 


356(7) 


0.624(73) 


4.912(89) 


1.027(13) 


0.546(14) 


M3 


356(7) 


0.648(67) 


4.944(97) 


1.051(13) 


0.557(14) 


M4 


495(10) 


0.787(25) 


3.908(35) 


1.056(7) 


0.566(7) 


M5 


597(12) 


0.814(13) 


3.495(17) 


1.074(4) 


0.574(4) 



Table IV. Full set of results for the scalar and tensor charges. The isovector data are renormalized in the MS scheme at 
(j, — 2 GeV, whereas the u + d data contain only the connected contribution and are not renormalized. Wilson data, W1-W8, 
are from the middle source-sink separation. 



where from the heavy quark limit 



| - Sa(l - m 2 <Sc L na), 



(B6) 



and the coefficient of the log term: 



6c 



LNA 



2(4^)2 



45l 
$9l 



NN 



NA \ 2 



9a 



(B7) 



Fixing mj, = 5 GeV as in |33j yields a fit with 3 independent parameters. In addition, fixing gA and F w as for the 
scalar charge, and ^""^ = 1.85, reduces the number of independent parameters to 2. 

Appendix C: Table of results 

Table ITvl contains the scalar and tensor data from our full set of ensembles. 
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